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The human LGN and SC consist of distinct layers, but their layer-specific response properties remain poorly 
understood. In this fMRI study, we characterized visual response properties of the magnocellular (M) and 
parvocellular (P) layers of the human LGN, as well as at different depths in the SC. Results show that fMRI is 
capable of resolving layer-specific signals from the LGN and SC. Compared to the P layers of the LGN, the M layers 
preferred higher temporal frequency, lower spatial frequency stimuli, and their responses saturated at lower 
contrast. Furthermore, the M layers are colorblind while the P layers showed robust response to both chromatic 
and achromatic stimuli. Visual responses in the SC were strongest in the superficial voxels, which showed similar 
spatiotemporal and contrast response properties as the M layers of the LGN, but were sensitive to color and 
responded strongly to isoluminant color stimulus. Thus, the non-invasive fMRI measures show that the M and 
P layers of human LGN have similar response properties as that observed in non-human primates and the super- 
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ficial layers of the human SC prefer transient inputs but are not colorblind. 


© 2015 Elsevier Inc. All rights reserved. 


Introduction 


A fundamental principle about the human visual system is the clear 
organization of anatomically and functionally distinct pathways or 
channels. In the geniculo-striate visual system, the two major pathways 
are the magnocellular (M) and parvocellular (P) pathways (Felleman 
and Van Essen, 1991), responsible for processing different aspects of 
the visual world (Kulikowski and Tolhurst, 1973; Legge, 1978). The M 
and P pathways are best represented and clearly segregated in different 
layers in the LGN. The inner two layers of the LGN consist of larger cells 
and are called the M layers, while the outer four layers with smaller 
cell bodies are the P layers. Although the physiological characteristics 
of M and P neurons in the LGN have been extensively studied in non- 
human primates (Derrington and Lennie, 1984; Jones, 2007; Wiesel 
and Hubel, 1966), little is known about the functional properties of 
the M and P neurons in humans due to the lack of non-invasive methods 
applicable to human LGNs. 
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The superior colliculus (SC) is another subcortical relay in the retino- 
tectal visual pathway. It is evolutionarily more primitive than the LGN 
and supports different but critical visual functions, such as head and 
eye movement, visual attention and even blindsight (Ro et al., 2004; 
Stoerig and Cowey, 1997). In primates, the SC is anatomically and func- 
tionally divided into superficial, intermediate, and deep layers. The 
intermediate and deep layers of SC are closely related to the controls 
of head and eye movement, receiving cortical inputs from a number of 
frontal and parietal cortical areas. The superficial layer of SC strongly 
responds to visual stimuli and receives direct input from the retina, in 
particular of the parasol ganglion cell type (Crook et al., 2008). The 
physiology of SC neurons has been extensively studied in nonhuman 
primates. In humans, the laminar response properties of the SC to visual 
stimulation remain to be explored. 

With the development of functional magnetic resonance imaging 
(fMRI) techniques, it becomes possible to measure the layer-specific re- 
sponse of human LGN and SC (Schneider, 2011; Schneider et al., 2004; 
Zhang et al., 2014), but it remains challenging due to their small sizes 
and the poor signal-to-noise ratio deep in the brain. In the first part of 
this study, we were able to reliably measure signals from the M and P 
layers of the human LGN. The identified spatial patterns of the M and 
P layers showed high consistency with the LGN anatomy. The functional 
properties of the human LGN measured with fMRI matched the evi- 
dence from monkey physiology studies (Derrington and Lennie, 1984; 
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Merigan and Maunsell, 1993; Wiesel and Hubel, 1966), with the M 
layers of the LGN tuned to higher temporal frequency, lower spatial 
frequency, and lower contrast, and vice versa for the parvocellular 
subdivision. While the P layers responded strongly to color stimulus, 
the M layers of the LGN were colorblind. In the second part of the 
study, we characterized the visual response properties of the human 
SC at different depth from the SC surface. Responses in the SC were 
strongest in the superficial voxels, which showed similar spatiotempo- 
ral and contrast responses as the M layers of the LGN, but were sensitive 
to color and strongly responded to isoluminant color stimulus. 


Materials and methods 
Subjects 


Thirteen healthy observers (age: 22.9 + 2.5 years, 4 males, 9 
females) participated in all the experiments, with normal or corrected 
to normal vision. All subjects were given written informed consents 
before the experiments. The institutional review board of Institute of 
Biophysics approved the experimental protocol. 


Stimuli 


Stimuli were generated in MATLAB (Mathworks Inc.) with psy- 
chophysics toolbox extension (Brainard, 1997; Pelli, 1997), running 
on a MacBook pro computer, presented by MRI safe projectors 
(1024 x 768 at 60 Hz). 

Visual stimuli were full-field (24° in diameter) counterphase flicker- 
ing checkerboard patterns, each presented for 16 s, alternating with 16 s 
of fixation. Stimuli were designed to preferentially engage the P and the 
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M systems, and they are labeled as the P stimulus and M stimulus for 
simplicity. As shown in Fig. 1A, the P stimulus was high spatial frequen- 
cy, isoluminant red/green square wave patterns, and reversed contrast 
at 0.5 Hz. The M stimulus was low spatial frequency sine wave patterns 
with 30% luminance contrast, and counterphase flickered at 10 Hz. For 
the P stimulus, the luminance of the red and green were matched for 
each observer with the minimal flicker procedure, the adjustment of 
red/green isoluminant was done 20 times, and the results were aver- 
aged. Fig. 1B shows the chromatic and achromatic stimuli. The chromat- 
ic stimulus has the same chromatic pattern as the P stimulus; the 
achromatic stimulus has the spatial frequency and was presented at 
50% luminance contrast. Both chromatic and achromatic stimuli 
counterphase flickered at 5 Hz. Fig. 1C shows the stimuli measuring 
the spatiotemporal frequency response functions. The luminance de- 
fined checkerboard patterns were manipulated at three different spatial 
frequencies (low, medium, and high SF), each presented at five different 
temporal frequencies (0.9375, 3.75, 7.5, 10, and 15 Hz), with contrast 
fixed at 50%. Fig. 1D shows the stimuli measuring the contrast response 
functions. The luminance contrast was manipulated at 4 different levels 
(6.25%, 12.5%, 25%, and 50%), and counterphase flickered at 7.5 Hz. All 
stimuli were presented at a mean luminance of 84 cd/m?, rotating clock- 
wise or counterclockwise at 2.4°/s (except for the M stimulus, which ro- 
tated at 24°/s) to prevent the signal drop from local adaptation. During 
the experiments, subjects were instructed to keep fixation and passively 
viewed the stimuli. 

All subjects participated in three scanning sessions performed in 
3 days. Each session consists of multiple runs. In the first session, 4 
runs were scanned for the M/P stimuli, and 4 runs for the chromatic/ 
achromatic stimuli, each run consists of 130 volumes of images. In the 
second scanning session, nine runs were scanned for stimuli with 
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Fig. 1. Visual stimuli. (A) The M and P stimuli. (B) The chromatic and achromatic stimuli. (C) The stimuli measuring the spatiotemporal frequency response functions. (D) The stimuli mea- 


suring the contrast response function. 
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different spatiotemporal frequencies, with 162 volumes for each run. 
Each run consists of all five temporal frequency conditions, and three 
runs were scanned for each spatial frequency condition. The third ses- 
sion measured the contrast response function, eight runs were scanned 
and 130 volumes were acquired for each run. 


MRI protocol 


MRI data were acquired with a 3 T MRI scanner (Siemens Trio) at the 
Institute of Biophysics, using a 12-element Head Matrix coil. A gradient 
echo planar imaging sequence was used to acquire functional images 
(2 mm isotropic voxels: 26 axial slices, 2 mm thickness, 128 mm FOV, 
64 x 64 matrix. TR/TE = 2000/28 ms, flip angle = 90°). High- 
resolution anatomic volume was obtained with aT1 MPRAGE sequence 
(1 mm isotropic voxels, 192 sagittal 1 mm thick slices, 256 x 256 matrix 
with 1 mm in-plane resolution, TR/TE = 2600/3.02 ms, flip angle = 8°). 
Another high-resolution T2 anatomical scan coaligned with the func- 
tional images was acquired to assist the registration of the EPI functional 
images to the T1 anatomical volume (TSE T2w sequence, TR/TE = 3400/ 
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Axial 


89 ms, 26 axial slices with 2 mm thickness, 256 x 256 matrix with 1 mm 
in-plane resolution, flip angle = 60°). 


Data analysis 


MRI data were analyzed using AFNI (Cox, 1996), SPM (Friston, 
1996), and custom MATLAB code. The first two volumes of each run 
were discarded to avoid magnetic saturation effects. The functional im- 
ages from each run were first resampled to 1 mm isotropic, realigned to 
the first image of the first run in this session. The high-resolution T1 vol- 
ume was registered to the T2 anatomical volume and then normalized 
to the MNI space. Functional images were normalized to the MNI 
space based on the warping parameters estimated from the normaliza- 
tion of the T1 volume. Linear trend in the fMRI time series were re- 
moved, and the unit of the time series was converted to percent signal 
change. We did not do any correction for physiological noise. According 
to previous studies, such correction should improve the SNR from 
subcortical nucleus (Limbrick-Oldfield et al., 2012). 


Sagittal 


Fig. 2. (A) Location and shape of the human LGNs on brain slices obtained from autopsy. Red arrows indicate the locations of the LGN. (B) Upper row clearly shows the shape of LGNs from 
T1 images of a representative subjects. Lower panels show the ROIs of the LGN for analysis. (C) The locations of the LGN in the full brain. Green crosses from different views indicate the 
same spatial location. (D) The LGNs from proton density weighted images. The LGNs can be clearly identified but individual layers are not resolvable. (E) ROIs of the SC were defined from 


T1 anatomical images. Here shows the coronal slices of a representative subject. 
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The ROIs of the LGN were defined based on combined functional 
activations from all stimuli of all sessions and the anatomical locations 
of the LGN from the T1 volume (second row in Fig. 2B). The LGNs can 
be clearly identified from the T1 mprage images. As a standard for 
comparison, Fig. 2A shows the human LGN on brain slices obtained 
from autopsy. On the T1 images in Fig. 2B, showing the LGNs from a rep- 
resentative subject, the LGNs appeared darker compared to surrounding 
brain tissues. The shapes of the LGNs from the T1 images (Fig. 2B, upper 
row) are highly consistent with those seen in the autopsy slices 
(Fig. 2A). Fig. 3C shows the location of the LGNs in the whole brain 
(Green crosses from different views correspond to the same spatial 
coordinate). In an attempt to resolve individual layers of the LGN, we 
collected high-resolution proton density images (512 by 512 matrix, 
0.5 mm in-plane resolution, with an 8-channel volume coil on a 7 T 
scanner) from one subject (Fig. 2D). The shape outline of the LGNs 
looks clearer than that in the T1 images, but individual layers are not 
resolvable. 

The beta values of different stimulus conditions were fitted and ex- 
tracted with general linear models (GLM). Beta values in the LGNs 
were extracted and analyzed in MATLAB. The coordinates of the right 
LGNs were mirror flipped to the left, and the mass centers of all LGNs 
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were registered. Group average of the beta values and functional T 
maps were generated from the registered LGNs. Similar approaches 
were used to analyze the data from the SC. 


Results 
Part 1: visual responses in the M and P layers of the human LGN 


We first identified the magnocellular and parvocellular layers of the 
LGN with stimuli designed to differentially activate the magnocellular 
(M stimulus) and parvocellular (P stimulus) neurons. Following the 
identification of these ROIs, we examined their responses to chromatic 
and achromatic stimuli, measured their spatiotemporal frequency 
response functions, and contrast response functions. 


Functional localization of the M/P layers of the human LGN 

As shown by Fig. 3A, the M and P stimuli were presented for 16 s, 
alternating with a 16 s period of fixation. Subjects were asked to keep 
the fixation during the experiment. For each LGN, the fMRI responses 
(beta values) to the M/P stimuli were fitted and extracted via general 
linear model (GLM). The LGNs were identified from significantly 
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Fig. 3. (A) Stimuli used for functional localization of the LGN's M and P layers. The P stimulus (0.5 Hz) and M Stimulus (10 Hz) were presented in separate 16 s blocks, alternating with a 16 s 
fixation period. (B) Topographies of the M and P layers of the LGN. The first column from the left illustrates the location of the LGN in the brain from coronal, axial, and sagittal views. The 
second column shows the enlarged spatial pattern of M/P layers identified by contrasting between fMRI responses to the M and P stimuli (paired t-test across all LGNs). The statistic map 
was paired t-test results across all LGNs (group results), over-layered on the anatomy of a representative subject. Purple lines indicate the boundary of the LGN volume. The color bar in- 
dicates the T values, only significant voxels were shown in the T map (Red: M beta > P beta, Blue: M beta < P beta). The third column shows the LGN anatomy as a reference, derived from 
human brain autopsy (the M layers were shown in red, and the P layers in blue color). The fourth column shows the locations of mass centers for the M/P functional subdivisions in coronal/ 
axial/sagittal planes. Error bars in this and all subsequent figures indicate standard errors of the mean. 
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activated voxels by the checkerboard stimuli at the anatomical locations 
of the LGN based on high-resolution T1 images. During group analysis, 
all LGNs (n = 26) were registered to their center of mass (right LGNs 
were mirror flipped to the left before registration), and a voxel-wise T 
test was performed to compared the difference between fMRI responses 
to the M and P stimuli. Fig. 3B (2nd column) shows the results. The M 
layers of the LGN were identified as voxels having significantly stronger 
responses to the M stimulus than to the P stimulus, and vice versa for 
the identification of the P layers of the LGN. Compared to the center of 
the P portion of the LGN, the center of the M portion was located signif- 
icantly more anteriorly (t(25) = 8.57, p < 0.001), inferiorly (t(25) = 
7.27, p < 0.001), and slightly more medially (t(22) = 1.48, p < 0.16) 
(Fig. 3B, 4th column). Individual coordinates for the M and P centers 
were shown with scatter plots in Supplemental Material. This spatial 
pattern is consistent with the anatomy of human LGN (Fig. 3B, 3rd 
column, derived from brain slices obtained by autopsy). 


Functional properties of the M/P layers of human LGN 

In separate fMRI scanning sessions, we measured the magnocellular 
and parvocellular neural responses of the human LGN along different 
stimulus dimensions. Fig. 4A shows the layer-specific fMRI responses 
of the LGN to the M and P stimuli. In the magnocellular layers of the 
LGN, the fMRI response to the M stimulus was significantly stronger 
than that to the P stimulus (t(25) = 4.97, p < 1074). While in the 
parvocellular layers, the fMRI responses showed a significant bias to 
the P stimulus (t(25) = 4.79, p < 1074). The fMRI responses of the M 
and P layers to chromatic and achromatic stimuli were shown by 
Fig. 4B. A significant two-way interaction was found between layers 


(M and P) and stimuli (chromatic/achromatic), F(1,25) = 32.6, 
p < 1077, indicating a dramatic difference in chromatic sensitivity of 
the magnocellular and parvocellular neurons. The M layers of the LGN 
showed strong fMRI response to the achromatic contrast stimulus but 
minimal response to the isoluminant chromatic contrast stimulus, 
whereas the P layers showed robust fMRI responses to both achromatic 
and chromatic patterns. Fig. 4C shows the temporal frequency response 
functions of the M/P layers. Compared to the P layers, the M layers 
showed stronger fMRI response to high temporal frequencies (10 Hz, 
15 Hz) but weaker response to low temporal frequencies (.94 Hz, 
3.75 Hz). Also, Fig. 4D shows that, compared to the P layers, the M layers 
had stronger response to the low and medium spatial frequency 
patterns but weaker response to high spatial frequency patterns. 
These results are consistent with the spatiotemporal response proper- 
ties of the M/P neurons. Fig. 4E shows the luminance contrast response 
functions. At the lowest contrast level (6.25%), the fMRI response of the 
M layers almost saturated and remained flat until 50% of luminance 
contrast. The P layers had weaker response than the M layers at the 
lowest contrast but showed more linear modulation of fMRI responses 
across different contrast levels. Our data are consistent with the LGN 
physiology that the M cells are more sensitive to low luminance contrast 
than the P cells. 

Thus, using fMRI, we were able to identify the magnocellular and 
parvocellular layers of the human LGN with stimuli preferentially acti- 
vating the magnocellular and parvocellular neurons. The spatial pattern 
of the identified M/P layers agrees well with the known anatomy of the 
human LGN, and the response properties of the M/P layers are also 
consistent with previous physiological evidence. 
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Fig. 4. Functional properties of the M and P layers of the human LGN. (A) fMRI response of the M/P layers of the LGN to the M/P stimuli. Red and blue bars indicate BOLD signals to the M and 
to the P stimulus, respectively. (B) The fMRI response of the M/P layers to the chromatic (red bars) and achromatic (blue bars) stimuli. (C) The fMRI responses as a function of temporal 
frequencies of flickering checkerboard. Red and blue lines show the fMRI responses in the M and P layers of the LGN, respectively. (D) Spatial frequency response functions of the M/P 
layers. (E) Contrast response functions of the M/P layers. The mean BOLD responses were normalized across subjects: subtracting the mean for each subject and then adding the 
grand-mean (across all subjects) back. Such normalization does not affect within-subject statistics, and the error bars better characterized the reliability of within-subject effects. The 


same normalization was done to generate Fig. 6. 
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Part 2: visual responses at different depths in the human SC 


We measured fMRI responses of the SC at different depths of the nu- 
cleus to a variety of visual stimulations, such as checkerboards defined 
by chromatic or luminance contrast, presented at different spatiotem- 
poral frequencies, and varied luminance contrast. 

The ROIs of the SC were defined for each individual according to their 
T1 anatomical images. The fMRI signals in the SC were averaged across 
all stimulus conditions, and the beta values to the averaged response 
were fitted with GLM. One sample t-test was performed on the beta 
maps across all SCs. Fig. 5A shows the group averaged beta map in the 
SC at significant voxels (t(25) = 4.0, p < 0.0002). The strongest response 
in the SC was located more laterally and superficially. The surface of the 
SC was defined as the outer most voxels, according to the anatomy of the 
SC. Then the depth for each voxel was calculated as the shortest distance 
from the voxel to all surface voxels. The third column in Fig. 5A shows 
the map of depth from the SC surface. As shown by Fig. 5B, fMRI re- 
sponse of the SC to the visual stimuli was strongest at the superficial 
voxels and gradually decreased with increasing depth of voxels from 
the surface. This result is consistent with the fact that the superficial 
layer of the SC is mostly composed of visual sensory neurons. 

Fig. 6 shows the functional responses of the SC voxels at different 
depth from the surface along multiple stimulus dimensions. As shown 
by Fig. 6A, voxels at the surface of the SC (depth = 0 mm) showed 
significant stronger responses to the M stimulus compared to the P 
stimulus (t(25) = 2.30, p < 0.03), and such difference was not signifi- 
cant deeper in the SC. It is unlikely that the signals from the superficial 
layers were from large veins, as the closest large blood vessel is located 
at least 3 mm away (to the medial side of the SC) from the center of SC 
activation (more lateral). Fig. 6B shows that voxels at the SC surface and 
those close to the surface responded significantly stronger to the achro- 
matic contrast stimulus compared to the chromatic contrast stimulus 
(depth = 0: t(25) = 2.30, p < 0.03, and depth = 1: t(25) = 2.30, 
p < 0.03). Importantly, these results showed that the SC strongly 
responded to isoluminant stimulus defined by chromatic contrast, 
suggesting that the SC is not colorblind. Figs. 6C and D show that the 
SC responses are better tuned to high temporal frequency and low 
spatial frequency stimuli. The contrast response functions of the SC 
were shown by Fig. 6E, the SC responses saturated at the lowest contrast 
level and remained almost constant across higher contrast levels. These 
results indicate that the spatiotemporal response properties, and lumi- 
nance contrast sensitivities of the SC are very similar to the functional 
properties of the magnocellular layers of the LGN. 
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To summarize, our data show that SC response to the visual stimuli 
was strongest in the superficial voxels and gradually decreased in 
deeper voxels. The human SCs responded robustly to stimuli defined 
by chromatic contrast, although luminance defined stimuli induced 
slightly stronger response than chromatic stimuli. Finally, human SCs 
prefer high temporal frequency and low spatial frequency stimuli, and 
their responses saturated at low luminance contrast. 


Discussion 


Although fMRI studies have made extensive progress at delineating 
the two-dimensional “functional areas,” relatively little attention has 
been directed to the third dimension—layers. The ability to access 
layer-specific neural signal is critical to the understanding of neural 
mechanisms in the brain. For example, a recent fMRI study found lami- 
nar differences in the neurovascular coupling in primate V1, providing 
important implications for the mechanism of feedforward and feedback 
processes (Goense et al., 2012). In this fMRI study, we were able to mea- 
sure BOLD signals from magnocellular and parvocellular layers of the 
human LGN as well as from different depth in the human SC. Their 
layer-specific visual response properties, such as chromatic selectivity, 
spatiotemporal response functions, and contrast response functions, 
were characterized. 

With stimuli designed to differentially activate the magnocellular or 
parvocellular neurons, the fMRI data reliably localized the M and P 
layers of the human LGN. The mass center for the M layers of the LGN 
were located more anteriorly, inferiorly and medially compared to the 
mass center for the P layers, which is consistent with the laminar profile 
of the human LGN. In addition, the fMRI response in the M layers 
favored higher temporal frequency, lower spatial frequency stimuli, 
and saturated at lower contrast, compared to the response in the P 
layers. Further, the M layers of the LGN did not respond to stimuli 
defined purely by chromatic contrast, while the P layers strongly 
responded to both chromatic and achromatic stimuli. These functional 
characteristics are in good agreement with the findings from previous 
monkey physiology studies (Derrington and Lennie, 1984; Merigan 
and Maunsell, 1993; Wiesel and Hubel, 1966). 

One previous study tried to identify the M/P layers of human LGN 
based on different contrast sensitivities of the M/P neurons (Schneider 
et al., 2004), but the segmented layers based on a single functional cri- 
terion showed large variability across subjects. Another study varied 
multiple feature dimensions in the localizer stimuli to differentiate the 
responses of the M/P neurons, and the identified the M/P layers did 
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Fig. 5. (A) Beta maps in the SC to visual stimulation. The first column shows the anatomical location of the SC in coronal, axial, and sagittal sections of the brain. The second column shows 
the group averaged beta map in the SC (map threshold is p < 0.0002). The third column shows the map of depth from the SC surface. (B) Visual responses in the SC as a function of depth 


from the surface of the SC. Signals were averaged across all stimulus conditions. 
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Fig. 6. Visual response properties at different depth in the human SC. (A) SC responses to the M/P stimuli. Red and blue bars indicate BOLD signals to the M and to the P stimulus, respec- 
tively. (B) SC responses to the chromatic (red bars) and achromatic (blue bars) stimuli. (C) SC responses as a function of temporal frequencies of flickering checkerboard. Colored lines 
indicate fMRI signals from different depth in the SC. (D) Spatial frequency response functions at different depth. E. Contrast response functions at different depth. 


show better consistency (Denison et al., 2014). However, their results 
came from a relatively small group of subjects (n = 6), and there 
were no independent experiments to investigate the functional proper- 
ties of the resultant M/P layers. In our study with 13 subjects, the M and 
P stimuli were optimally designed along four stimulus dimensions (con- 
trast, color, spatial, and temporal frequencies) to selectively favor the 
responses of the M or P layers of the LGN. The topography of the identi- 
fied M and P layers was in good agreement with the anatomy of the 
human LGN. The robustness of this approach in identifying LGN layers 
was supported by the observation that essentially the same M and P 
topography was obtained from a separate group of 18 subjects scanned 
at a different scanner with the same stimuli and procedure (Zhang et al., 
2014). Furthermore, in a series of independent fMRI experiments, we 
characterized the functional properties of the M and P layers, which 
are consistent with the LGN physiology. Therefore, the current study 
provided convincing segmentation results of the M/P layers of human 
LGN, as well as detailed characterization of their visual response 
properties. 

A more conclusive approach might be to identify the boundary be- 
tween the M and P layers of the LGN with ultra-high resolution MRI 
techniques such as proton density imaging sequences (McKetton and 
Schneider, 2012), and then validate whether the M and P functional 
subdivisions conform to that boundary. We tried to collect high- 
resolution PD images of the LGN from one subject. As shown by Fig. 2, 
the LGN body can be better visualized than the T1 images, but individual 
layers remain unresolved. It is possible that using a bite bar during scan- 
ning may help. However, in the current study, the LGN bodies are clearly 
visible from the T1 anatomical images, and the ROI of the LGN were 
selected as voxels within the LGN body that were significantly activated 
by the visual stimulation. Therefore, the identified M and P functional 


sections are within the anatomical volume of the LGN. In addition, con- 
sidering that the relative locations of the M and P sections are consistent 
with the known anatomical relationships, we have high confidence 
about the current segmentation results. 

In the second part of this study, we characterized visual response 
properties at different depths in the human SC. Results show that visual 
responses in the SC were strongest in the superficial voxels and gradu- 
ally decreased deeper into the nucleus. This finding is consistent with 
physiological results obtained in non-human primates that the superfi- 
cial layer of the primate SC is mostly composed of visual sensory 
neurons. Both chromatic and achromatic stimuli induced significant 
fMRI response in the SC, suggesting that the human SC is not colorblind. 
A recent monkey single-unit study found that chromatic stimuli evoked 
spiking responses of neurons in the intermediate layers of the SC (White 
et al., 2009), whereas our human fMRI data showed strongest fMRI 
responses to chromatic stimuli in the superficial layers of the SC. Such 
discrepancy could be explained by the different neural basis for the 
fMRI signal and single-unit spiking activities. Compared to neuronal 
spiking activities, the fMRI BOLD signal might be more sensitive to 
top-down feedback modulations. In fact, previous studies showed that 
attention enhancement for visual responses (in contrast to visuomotor 
activities) in the human SC are strongest in the superficial voxels 
(Katyal and Ress, 2014; Katyal et al., 2010). Finally, fMRI response in 
the SC preferred high temporal frequency and low spatial frequency 
stimuli and saturated at low luminance contrast. These visual response 
properties are similar to what we found for the M layers of the LGN. To 
summarize, our results are consistent with previous fMRI studies on the 
human SC that showed strongest visual responses in the superficial 
layers (Katyal et al., 2010) and low sensitivity to luminance contrast 
(Schneider and Kastner, 2005). Importantly, our data also revealed the 
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spatiotemporal response properties of the human SC, which was similar 
to responses in the M layers of the LGN, and showed a counterintuitive 
finding that the human SCs are very sensitive to color. 

The subcortical pathways at the early stages of visual information 
processing have far reaching consequences and broad implications. 
The ability to access the layer-specific responses of the human LGN 
and SC is of great value to study the broad roles of these early visual 
pathways in human visual perception and cognition. For example, it 
could allow us to address issues such as whether visual consciousness 
is more related to one or the other pathway. It can also help to investi- 
gate important clinical issues, including the differential involvement of 
the parallel pathways in certain neurodegenerative diseases such as 
glaucoma, developmental dyslexia, and amblyopia. 

Supplementary data to this article can be found online at http://dx. 
doi.org/10.1016/j.neuroimage.2015.02.025. 
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